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Abstract. Nonpolynomial basic hypergeometric eigenfunctions of the Askey- Wilson second 
order difference operator are known to be expressible as very-well-poised s4>7 series. In this 
paper we use this fact to derive various basic hypergeometric and theta function identities. 
We relate most of them to identities from the existing literature on basic hypergeometric 
series. This leads for example to a new derivation of a known quadratic transformation 
formula for very-well-poised g07 series. We also provide a link to Chalykh's theory on (rank 
one, BC type) Baker- Akhiezer functions. 
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1 Introduction 

In this paper we present derivations of various basic hypergeometric and theta function iden- 
tities using the interpretation of very- well-poised 8<^7 series as eigenfunctions of the Askey- 
Wilson second order difference operator T>. For instance, we reobtain a nonstandard type three 
term transformation formula for very-well-poised 8^7 series [7] as the connection formula for the 
asymptotically free eigenfunctions of T>, we investigate the eigenfunctions of T> with trivial quan- 
tum monodromy and relate them to (rank one, BC type) Baker-Akhiezer functions from [4, 5], 
we rederive the quadratic transformation formula [6, (3.5.10)] for very-well-poised 8</>7 series, 
and we obtain various theta function identities by translating symmetries of the Askey-Wilson 
function and of the asymptotically free eigenfunction of T> in terms of the associated normalized 
c-functions. 

The most general family of orthogonal polynomials satisfying a second order g-difference 
equation is the family of Askey-Wilson polynomials [1]. For our purposes it is convenient to 
view the associated difference equations as eigenvalue equations for the second order Askey- 
Wilson second order difference operator T> already mentioned in the previous paragraph. The 
operator V depends, besides on the difference step-size and the deformation parameter q, on 
four additional free parameters. 

Nonpolynomial basic hypergeometric eigenfunctions of T> have been subject of study in va- 
rious papers (see, e.g., [8, 9, 10, 19, 20, 24, 26]). Important examples are the Askey-Wilson 
function £(-,z) and the asymptotically free eigenfunction <&(-,z) (the corresponding eigenvalue 
depends in an explicit way on q z + q~ z G C). They are defined provided that < \q\ < 1. 
They are selfdual eigenfunctions (the role of the argument and z is interchangeable). They 
naturally arise in harmonic analysis on the quantum SU(1, 1) group and in the study of the 
double affine Hecke algebra of type C v Ci (the rank one Koornwinder case), see, e.g., [11] 
and [24, 25] respectively. From this representation theoretic viewpoint £(-,z) plays the role of 
the spherical function and $(•, z) the role of the Harish-Chandra series. 
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Ruijsenaars' i?-function [19] is another nonpolynomial selfdual eigenfunction of T> which is 
required to satisfy yet another second order difference equation of Askey-Wilson type. The 
step-direction of the two Askey-Wilson second order difference equations are allowed to be co- 
linear (which corresponds to deformation parameter q being of modulus one). The R- function 
arises as matrix coefficient of representations of the quantum double of the quantized universal 
enveloping algebra of sfe, see [2]. 

The Askey-Wilson function £ (•, z) can be explicitly expressed in terms of the asymptotically 
free eigenfunctions <J>(-,±£). The elliptic function c(-,z) governing the expansion coefficients, 

£{; Z) = C(; *)*(■, z) + C(; -«)*(-, -z), 

is called the (normalized) c-function [10]. It is explicitly given as quotient of theta functions. 
The selfduality of £ and $ and the fact that £(—x, z) = £(x, z) then allow us to express x, z) 
in terms of $(a;, ±z) (connection formula). The cases when the connection coefficient formula 
trivializes is particularly interesting since it directly relates to the theory of Baker-Akhiezer 
functions [4, 5]. We discuss this in Section 3. 

Suitable two parameter specializations of the Askey-Wilson polynomials yield the continuous 
g-Jacobi polynomials. They have appeared in two guises (see [17] and [1]) which are interrelated 
by a quadratic transformation formula going back to Singh [22], see also [1, § 4] and [6, § 3.10]. 
This quadratic transformation formula was derived in [1] using the orthogonality relations of the 
Askey-Wilson polynomials. Ruijsenaars [20] stressed that for these parameter specializations 
the Askey-Wilson second order difference operator T> factorizes up to an additive constant as 
a square of an Askey-Wilson type second order difference operator with step-size half of the step- 
size of T>. We use this observation to prove a quadratic transformation formula for a two param- 
eter family of the asymptotically free eigenfunction $(•, z) of T>. This complements Ruijsenaars' 
results [20], where he lifted the quadratic transformation formula for continuous q-Jacobi poly- 
nomials to a quadratic transformation formula for a two parameter subfamily of the i?-function. 

Using the known explicit expression of §>(-,z) in terms of basic hypergeometric series we 
link the above mentioned results to various known identities for very-well-poised s(f>7 series. For 
example, the quadratic transformation formula for <£(•, z) becomes the known quadratic transfor- 
mation formula [6, (3.5.10)] for very- well-poised 8</>7 series after applying suitable transformation 
formulas to both sides of the identity, see Remark 5.3(a) (this was observed by M. Rahman). 

Combining symmetries of the Askey-Wilson function with its c-function expansion yields 
nontrivial identities for the c-function, hence nontrivial theta function identities. For instance, 
the fact that £(-,z) is invariant under negating the argument yields a theta function identity 
which is a six parameter dependent subcase of a theta function identity [6, Exercise 5.22] due to 
Slater [23]. The fact that $(•, ±z) satisfies a quadratic transformation formula but the Askey- 
Wilson function £ (•, z) does not, yields a nontrivial identity for the c-function and consequently 
a quadratic type theta function identity (7.4). 

The results in the present paper play an important role in the study of the spectral problem 
of the trigonometric Macdonald-Ruijsenaars-Cherednik commuting family of difference ope- 
rators associated to root systems. They are for instance needed in the asymptotic analysis 
of g-analogues of Harish-Chandra series associated to root systems (cf. [13, 14, 15, 25]) along 
codimension one facets of the Weyl chamber. I will return to this topic in a future work. 

2 Eigenfunctions of the Askey-Wilson 
second order difference operator 

We use, besides a deformation parameter < q = e 27rv/ ~* T < 1 (r € \J— lK>o) and a choice 
of step-size s G Q>0; four free parameters {k, X, v, which we call Hecke parameters (this 
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name comes from their interpretation as multiplicity parameters in the Cherednik-Macdonald 
theory of the double affine Hecke algebra of type C v Ci, see [16]). We assume that the Hecke 
parameters k, X, v and ? are real. We set q x := e 2 ~ K v r - : ^ TX for x £ C. In the paper q will be 
fixed throughout. Step-sizes s and | will simultaneously appear in formulas. One can restrict 
without loss of generality to considering step-sizes 2 and 1, but formulas are more transparent 
when arbitrary values s of the step size are taken into account because it makes the s-dependence 
of the Askey- Wilson parameters 

{a s ,b s ,c s ,d s } := {q K +\-q K -\qi+^,-qi+ v -"} 

explicit. We define dual Askey- Wilson parameters by 

{a s ,b B ,Cs,ds} ■= {q K+v ,-q K - v ,q^ +X+ ',-q^ +X -'}, 

i.e. the roles of the Hecke parameters A and v are interchanged. The parameters a s , b s , a s 
and b s do not depend on s, we therefore occasionally omit the subindex s for these Askey- 
Wilson parameters. Interchanging A and v defines an involution on the set of Hecke parameters, 
hence also on the associated set of Askey-Wilson parameters. In addition we have 

a 2 s = q~ s abcd, a s b s = a s b s , = c s d s , 

b s 

~~ _ Q s a s _ i ^ ~7- a g s "s _ u 

c s d s 

We now first recall the asymptotically free eigenfunction of the Askey-Wilson [1] second order 
g s -difference operator, which we will regard here as a second order difference operator with step 
size s. Explicitly, the Askey-Wilson second order difference operator T>, acting on meromorphic 
functions on C, is defined by 

(Vf) (x) := A(x) {f{x + s)- f{x)) + A{-x) (f(x - s) - f(x)) , 
;i - a s q x )(l - b s <f)(l - c s q x )(l - d s q x ) 



A{x) :-- 



a s (l - q 2x )(l - q s + 2x ) 



Sometimes it is important to write explicitly the dependence on the parameters, in which case 
we write T> as T>i s \ „, ,..„. 

Remark 2.1. The Askey-Wilson second order difference operator T> can be interpreted as 
a second order g s -difference operator when acting on r~ ^translation invariant meromorphic 
functions on C (which are the meromorphic functions of the form g(q x ) with g meromorphic 
on C*). 

We now define the elementary function W(x, z) = W(x, z; k, A, v, q, s) by 
W(x,z) := g(«+A+*)(«-HH-*)/«_ 

Note that W(x + s,z) = q K+v+z W{x, z) = aq z W(x, z). 

For generic bj the r +i</v basic hypergeometric series is defined by the convergent power series 



r+l9r 



ai ,a 2 ,...,a r+1 A » (ai,q a ,...,qr + i;g) j 
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s j-i 

where (ai, . . . , a s ; q) . = n Yl (1 ~~ Or<f) for j G Z>o U {00} (empty products are equal to one 



3 r =\ i=o 



by convention). The very- well-poised s</> 7 series is defined by 



1 1 

5 M/ 7 (ao;ai,a2,a3,a4,a5;Q , ,^J = 807 I 1 _i ; g, z 

\ao , -«o 2 3 90o/aii • • • ,gao/«5 / 

l-a q 2r ^ r yr [Uj;q) r 
I- an Z M 



r=0 



In case 2 = o;o { Z 2 / a i Q; 2«3tt4a5 it has a meromorphic continuation to (C*) 6 as function of 
{an 1 on, ■ ■ • ,05)- This follows from the identity [6, (III. 36)] expressing such very-well-poised ^7 
series as a sum of two 4^3 series. Define the holomorphic function St (a;) = St(x; k, A, v, q, s) 
("St" is standing for singular term) in x G C by 

St(x) := [q^/a a ,^/b.,q*yc.,q^/d s -,q>) oB . 
We write St d {z) := St{z; k, v, A, q; q, s) for the singular term with respect to dual parameters, 

St d {z) = {q s+z /a s ,q s+z /b s ,q s+z /c s ,q s+z /d.)?)^ 
The following proposition combines and refines observations from [9, 10, 25]. 

Proposition 2.2. There exist unique holomorphic functions T r on C (r > 0) satisfying the 
following three conditions, 

(1) T {z)={q^ 2z ;q s ) oo . 

00 

(2) The power series ^f{x,z) := ^ T r {z)q rx is normally convergent on compacta of {x,z) £ 

r=0 

C x C {consequently ^f{x, z) is a holomorphic function in {x, z) £ C x C). 

(3) The meromorphic function 



W{x,z) 
St{x)St d {z) 



Hx,z) := — \, *{x,z) (2.1) 



satisfies 

{ {V x - q z - q~ z + a + a' 1 ) $) {x, z) = 0, (2.2) 

(s) 

where T> x stands for the Askey-Wilson second order difference operator D = T> KXv(; . q 
acting on the x-variable. 

Furthermore, T r is r _1 -translation invariant and 

s+*+z as gs+x+z^ q s+x+*c s g s +"+ 2 a s s+ 2z J n x. n s\ 



n s+x+2z n s+z n s+z ____ „s+a; 

d s a s d s d s 



if\d s q x \ < 1. 
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Proof. The explicit expression (2.3) is 

( q s +*+*a s q a + x + z b s q a + x +*c s q s + x + z H s *+2z J n x. n s\ 
\ as ' a s ' a 3 ' d s >9 i «*« >« 



z) 



r=0 



r+2sr+x+2z x (g_ g^^^g—;^)^ 

do J ( 9 S+X+Z aa q s +x+ z as g s +*+*b a q=+^+-c s nS +2z. n s\ ' ^ ' ' 

\ d s ' a s ' a s ' a s 



It is a well defined meromorphic function in (x, z) G C x C provided that |d s g x | < 1, with 
possible poles at q s + sr + x + 2z = d s ( r g Z>o). It can be expressed as sum of two 4(^3 series using 
[6, (III. 36)] with parameters (a, 6, c, d, e, f, q) in [6, (III. 36)] specialized to 

(q s+x+2z /d s , q s+z /a s , q s+z /d s ,b s q Z ,c s q z , q s+x /d s , q s ), 

leading to the expression 

( g s+x j n x gflf sf^f. g 3+3:+z a s gS +x+ ZSa _ , 

/ f a s q x , q ^,b s q z ,c s q z s A 

X 493 I q °+x+Zg s q °+x+Zq s q"g a > 1 > 3 I 

\ o» ' ds ' d« / 

■f +x 1 z ~ z q s+x+z a s q s + x + z d a 



1-2 „S + Z 



x ^( ^4^£ ^^^)- (2 - 5) 

This alternative expression provides the meromorphic continuation and shows that *$>(x,z) is 
holomorphic in (x,z) € C x C. The expression (2.4) of ^(x,z) shows that $?(x,z) satisfies (1) 
and (2). By [9] (see also [10] for notations close to the present one), the resulting meromorphic 
function $>(x,z) (see (2.1)) indeed satisfies the difference equation (2.2). 
It remains to prove uniqueness. If a series of the form 

is a formal solution of (2.2) then the T r (z) (r > 0) satisfy recursion relations of the form 

r-1 

a( q sr _ 1){q z _ q sr-z )Tr{z) =J2v r t (q*)r t (z), r > 1 

t=0 

for some Laurent polynomials v\. This shows that the T r (z) (r > 1) are uniquely determined 
I'.v I oi :). ■ 

If confusion may arise about the parameter dependencies then we write the asymptotically 
free solution <£(x, z) as z\ k, A, v, <r; q, s) and *ff(x, z) as ^(x, z; k, A, v, q q, s). 

Remark 2.3. The characterization of 3>(x, z) as eigenfunction of T> is equivalent to a characteri- 

oo 

zation of ^f(x, z) = T r (z)q rx as eigenfunction of the second order difference operator obtained 

r=0 

from V by gauging it with gauge factor W(-, z)/St(-). The gauged difference operator and the 
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relevant eigenvalue are symmetric under arbitrary permutations of the Askey-Wilson parame- 
ters a s , b s , c s , d s . Since the normalization Tq(z) = [q s+2z ] q s ) of ^(x, z) is independent of the 
Askey-Wilson parameters, it follows that ^f(x,z) is symmetric under arbitrary permutations of 
the Askey-Wilson parameters a s , b s , c s , d s . 

Let 9(u;q) := (u,q/u;q) be the modified Jacobi theta function and write 

r 

6(m, . . .,u r ; q) = Jj0(iti;g) 

i=i 

for products of theta functions. Define the (normalized) c-function c(x, z) = c(x, z; k, A, v, q, s) 

by 

9(a s q- z ,b s q~ z ,c s q~ z , ck ^;<f) 

W(x,z)9(q Zz ,d s q x ;q s ) 

Using 9(qu; q) = —u~ 1 9(u; q) it follows that c(x + s, z) = c(x, z) and c(x, z + s) = c(x, z). We 
write c d (x, z) = c(x, z\ K, v, A, q q, s) for the c-function with respect to dual parameters. 
The Askey-Wilson function £(x, z) = £( A, v, q, s) is defined by 

£/ x z \ . = {~d a ' °^ds ,a s b s ,a s c s , 1 ^ L ;q s ) O0 



,^,± 7r ,a s b s c s q z ;q s ) c 



x 8 W 7 (^ s b s c s q' s+z ;a s q x ,a s q- x ,a s q z ,b s q z ,c s q z ;q s , ^ 

for \q s ~ z /d s \ < 1 (see [10]). Using [6, (III. 36)] to express £(x, z) as sum of two 4^3 series it follows 
that £(-,•) has a meromorphic extension to C x C. We write £ d and <3? d for the meromorphic 
functions £ and $ with respect to dual parameters. The following properties of £ and are 
known from [10, 25] (cf. also [8, 9, 26]). 

Proposition 2.4. 

(i) £(x, z) = £ d (z, x) (self duality) . 
(ii) £(—x,z)=£(x,z) and £(x,—z) = £(x, z). 
(Hi) 0>(x, z) = & d (z, x) (self duality) . 

(iv) £(x, z) = c(x, z)&(x, z) + c(x, —z)$(x, —z) (c-function expansion). 

Proof, (i) This follows from the transformation formula [6, (III. 23)] for very-well-poised 807 
series. 

(ii) By the explicit expression (2.7) it is clear that £(—x,z) = £(x,z). By (i) it then also 
follows that £ (x, —z) = £(x, z). 

(Hi) This follows again by application of the transformation formula [6, (III. 23)] for very- 
well-poised 8^7 series. 

(iv) Use Bailey's three term transformation formula [6, (III. 37)] for very-well-poised 8</>7 
series. ■ 

Remark 2.5. The selfduality of $>(x,z) and £(x, z) ensures that $>(x, ■) and £(x,-) are eigen- 
functions of the Askey-Wilson second order difference operator with respect to dual parameters. 

Note that the Askey-Wilson second order difference operator T> is invariant under x 1— > —x, 
hence 3>(— x, z) is again an eigenfunction of T> x with eigenvalue q z + q~ z + a + aT l . It can be 
expressed in terms of the eigenfunctions $>(x, z) and <3?(x, —z) as follows. 
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Corollary 2.6. 

fc(x, z) — c d (z, x) \ T . , c(x,—z) . . 

$(-x,z) = ( V J, Y-^- )$(x,z ) + -tt 2 2.8 

\ c d (z,-x) J c d (z,-x) 1 ' ; y ' 

(connection formula). 

Proof. By the previous proposition we have 

£{x, z) = £ d (z, -x) = c d (z, -x)<f> d (z, -x) + c d (z, x)¥(z, x) 
= c d (z, x)$(x, z) + c d (z, -x)$(-x, z). 

Compared with the c-function expansion for £(x, z) (see Proposition 2A(iv)) we get the desired 
result. ■ 

Remark 2.7. The connection formula (2.8) is not a direct consequence of Bailey's three term 
transformation formula [6, (III. 37)] for very-well-poised 8<^7 series. This is reflected by the fact 
that the coefficient of $>(x, z) in (2.8) does not admit an explicit expression as a single product 
of theta functions. The connection formula (2.8) is though directly related to the three term 
transformation formula [7, (5.8)]. 

Corollary 2.8. The c-function satisfies 

c(x, z)c d (z, -x)c d (-z, -x) + c(— x, z)c d (z, x)c d (-z, -x) 

= c(x, z)c(—x, z)c d (—z, —x) + c(x, z)c(—x, —z)c d (z, —x). (2.9) 

Proof. Since £(—x, z) = £(x, z) we have 

c(x, z)&(x, z) + c(x, — z)$(x, —z) = c(—x, z)&(—x, z) + c(—x, —z)&(—x, —z). 
Applying (2.8) twice to the right hand side of this formula implies 

a(x, z)$(x, z) = -a(x, -z)®(x, -z) (2.10) 

with the function a(x, z) given by 

c(— x, z)c(x, z) — c(— x, z)c d (z, x) c(—x,—z)c(x,z) 

a[X, Z) = T7 : 1 T, ; c(x, z). 

v ' ' c d (z,-x) c d (-z,-x) v ; 

Replace in (2.10) the variable x by x + ms and consider the asymptotic behaviour as m — > oo of 
both sides, using the fact that c(x, z) is s-translation invariant in both x and z and using that 

a-' m q- mz ^(x + ms, z) = T (z) (l + 0{q sm )) 

as m — > oo. It gives a(x, z) = 0. This is equivalent to (2.9). ■ 

Remark 2.9. Substituting in (2.9) the explicit expression (2.6) of the c-function c(x,z) gives 
the theta function identity 

tq x+z , tq x ~\ aq x ,bq X , cq x , dq- x , q 2z ; q) - 9 (iq~ x+z , tq x+z ,aq Z M ,cq z , dq~ z , q 2x ; q 

= q 2x 9 (iq~ x - z , iq~ x+z , aq~ x , bq~ x , cq~ x ,dq x , q 2z ; q 
\a a 

- q 2z 6 (tq- x - z , tq x - z ,aq- Z M-'\cq- z , dq z , q 2x ; q ) 
\ a a 
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where we have taken s = 1 and have written a = a\,...,d = d\ (and similarly for the dual 
parameters). This is a special case of Slater's theta function identity [6, Exercise 5.22], with the 
parameters (a,b,c,d,e,f,g,h) in [6, Exercise 5.22] specialized to 

( i, U-^M- z M-\~c q -\ C i<f+\ qA . 

V a d a J 

Note that Slater's formula [6, Exercise 5.22] is more general since it has, besides q, seven free 
parameters, while the formula (2.9) only has six. 

3 The case of trivial quantum monodromy 

Determining the monodromy representation for Gauss' second order hypergeometric differential 
equation is equivalent to deriving its connection coefficient formulas. The connection coefficient 
formulas explicitly relate the fundamental series expansion solutions around the three regular 
singularities of the hypergeometric differential equation. These formulas turn out to be di- 
rectly related to well known three term transformation formulas for the Gauss' hypergeomet- 
ric function %Fx. From the above notion of monodromy only its incarnation in terms of con- 
nection coefficient formulas generalizes to the difference setup (see Sauloy [21] and references 
therein for a detailed discussion of this issue). We therefore say that the explicit connection 
coefficient formula (2.8) solves the quantum monodromy problem of T>. In addition we say that 
the quantum monodromy is trivial if the coefficient of <&(x, z) in (2.8) vanishes. The latter 
terminology is motivated as follows. Firstly, for trivial quantum monodromy the connection 
coefficient formula (2.8) reduces to a simple equivariance property of the following renormaliza- 
tion <3>(x, z) = z; k, A, v, q q, s) of the asymptotically free eigenfunction <3?(x, z), 

<3?(x, z) := c(x, z)$(x, z), 

see Proposition 3.1(ra) below. Secondly, for an important subclass with trivial quantum mon- 
odromy, a suitable renormalization of the asymptotically free eigenfunction &(x, z) has a termi- 
nating series expansion, see Remark 3.2(z) for further discussions on this issue. 

Note that the renormalization $(x, z) of 3>(x, z) still satisfies the eigenvalue equation 

V(${-,z)) = (q z +q- z -a-a- l )l(;z) 

since c(x, z) is s-translation invariant in x. Write § d (z, x) := $(z, x; k, v, A, q, s). 
Proposition 3.1. 

(i) If j, j, | G with an even number of them being integers, then c(x, z) = c d (z, x). 

In the remaining items of the proposition we assume that (k, A, v, ?) is a four-tuple of real 
parameters such that c(x,z) = c d (z,x). 

{ii) $(x,z) = <& d (z, x) (selfduality) . 

(Hi) <&(—x,z) = <b(x,—z). 

(iv) £ (x, z) = <&(x, z) + $(— x, z). 

Proof, (i) Consider the quotient a := c(x, z)/c d (z, x) (we suppress the dependence on the 
variables and parameters). Using the explicit expression of the normalized c-function (2.6) in 
terms of theta functions it follows that a = lif(j,^, : |,|)is taken from the set 

{(0,0,0,0), (1, i,o,o),(i,o,i,o), (i,o,o,i), 

(°' \\ ■ °) - (°i °'£) - ("• 11) • {ii 11) } ■ 
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In addition, it is easy to check that a is invariant under integral shifts of the rescaled Hecke 
parameters j, |, |. This gives the result. 
(ii) This follows from Proposition 2.4(ra). 

(Hi) This follows from the connection formula (2.8) for <&(x, z). 

(iv) This is immediate from the c-function expansion of the Askey-Wilson function, see 
Proposition 2A(iv). M 



Remark 3.2. 

(i) Due to [13, § 4] the function $(x, z) for Hecke parameters (k, X,v, satisfying, besides 
the integrality conditions from Proposition 3.1(i), suitable additional (positivity) condi- 
tions, is up to normalization Chalykh's BCi type normalized Baker-Akhiezer function [4, 
§ 6]. An important property of the Baker-Akhiezer function is the fact that it has a ter- 
minating series expansion. Surprisingly this is not evident from the explicit expressions 
of <&(x,z) in terms of basic hypergeometric series. The properties from Proposition 3.1 
can be easily matched with the properties (see [4, 5]) of the Baker-Akhiezer functions. 
For instance, the Askey-Wilson function £(x,z) relates to the symmetrized normalized 
Baker-Akhiezer function <J>_|_ from [5, (3.20)]. In particular, formula [5, Theorem 3.9] re- 
lating the symmetrized Baker-Akhiezer function to Askey-Wilson polynomials matches 
with the known fact that the Askey-Wilson function £(x, z) reduces to the normalized 
Askey-Wilson polynomial for suitable discrete values of z, see (6.2) below. 

(ii) The relation between the Ai type Baker-Akhiezer function (see [4, § 4.1]) and Heine's 
basic hypergeometric series 2<f>i was stressed in an informal note of Koornwinder [12]. 



4 The factorization of the Askey-Wilson 
second order difference operator 

Ruijsenaars [20] analyzed when the square (T> + a + a -1 ) 2 is an Askey-Wilson second order 
difference operator again (with doubled step-size). An important special case turns out to be 
when the Hecke parameters are of the form (k, A, 0,0). In our notations the resulting formula 
[20, (3.11)] reads as 



[V 



(8/2) 

k,A,0,0;<? 



+ <f + q~ 



(T>(»/2) , 
V^K,A,0,0;g + 



+ q 2 +q K + q K ) 



V 



(s) 

K,X,K,\;q 



■2k 



+ ^ + q 



-2k 



(4.1) 



Remark 4.1. Our operator T> + a + a^ 1 = T> + q 2K + q~ 2fi for arbitrary Hecke parameters 
(k, A, v,s) is essentially the second order difference operator A$(c; •) in Ruijsenaars' paper [20, 
(2.1)] with the parameters (cq, ci, C2, C3) in [20] related to our Hecke parameters by k + A = 
— v— Icq, k — A = — V — lei, v + ? = — \/— Tc2 and v — ^ = —^J^lc^. 



A simple derivation of (4.1) is as follows. Consider the second order difference operator 
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Then 

Formula (4.1) thus is equivalent to 

(4:5 - 9 f - 9- f ) (4:5 + ? f + ) = <L,a ;? - - t 2 + ^ + ^- 2k > (4-2) 

which is an easy check (use that the Askey-Wilson parameters associated to Hecke parameters 
(re, A, re, A) satisfy c s = q^a s = gi+ K + A and d 8 = q^b s = —qz +K ~ x ). 

The eigenfunctions of T> x of the form p n (q x + q~ x ) (re G ^>o) with p n a polynomial of 
degree re, are the well known Askey-Wilson polynomials [1] (see also Section 6). In the special 
case that the associated four Hecke parameters are taken to be (re, A, re, A) (corresponding to the 
condition that the Askey-Wilson parameters satisfy c s = q^a s and d s = q^b s ) the Askey-Wilson 
polynomials are called the continuous g-Jacobi polynomials by Askey and Wilson [1, § 4]. In 
[1, (4.22)] a quadratic transformation formula for balanced 4^3 series, going back to Singh [22], 
is used to relate the above continuous g-Jacobi polynomials to Rahman's definition [17] of the 
continuous g-Jacobi polynomials (see also Section 6). 

Ruijsenaars [20] used the factorization (4.1) to motivate and analyze quadratic transforma- 
tion formulas for the hyperbolic nonpolynomial generalization of the continuous g-Jacobi poly- 
nomial, which is Ruijsenaars' i?-function with the continuous g-Jacobi specialization (re, A, re, A) 
of the associated Hecke algebra parameters (see Remark 4.1 for the relation with Ruijsenaars' 
notations [20]). In the following section we use the factorization (4.1) to prove quadratic trans- 
formation formulas for the asymptotically free eigenfunction <£(•, z; re, A, re, A; q, s). 



5 Quadratic transformation formulas 

For a function S(x, z; re, A, re, q q, s) we write 
E R (x, z) := E (x, |; re, A, 0, 0; q, 



This is the parameter specialization which reduces the Askey-Wilson polynomials to Rahman's 
version of the continuous g-Jacobi polynomials [17] (in base qz). Furthermore we set E R (x, z) := 
H(x, I ; re, 0, A, 0; q, |) for its dual version. Since 

we know from the previous section that the meromorphic function & R (x, z) satisfies 

{C x -q% -q-^ R (;z) =0. (5.1) 
In fact, <J?r(x, z) is the unique solution to (5.1) which is of the form 

R{ } St R (x)St d R (z) ^ R ' rK H 

with Y R j.{z) holomorphic in z G C, with T r< q(z) = (q^ +z ] q?) an d with the power series 

00 

ty R (x,z) := ]P T R)r (z)q rx converging normally on compacta of (x,z) G C x C. 

r=0 
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For a function H(x, z; k, A, v, q q, s) we write 

Ej(x, z) := E(x, z; k, A, «, A; g, s). 

This is the parameter specialization which reduces the Askey-Wilson polynomials to Askey's 
and Wilson's version of the continuous g-Jacobi polynomials [1, § 4]. In addition we write 
z) := H(x, z; k, k, A, A; q, s) for its dual version. 

Theorem 5.1. <&j(x,z) = $r(x,z). 

Proof. First of all note that the Askey-Wilson parameters associated to the Hecke parameters 
(re, A, 0, 0), deformation parameter q and step-size | are given by 

(q«+\- q «-\ q i,-qi). 
The expression of *$>(x, z) as a sum of two 4^3's, see (2.5), shows that 

(-qi+ x ;q i)^ R (x,z) 
is holomorphic in (x, z) € C x C. By Remark 2.3 also 

(ql+ x : -ql+ x ;ql)^y R (x,z) = (gf +2x ; q s )^V R {x, z) 
is holomorphic in (x, z) € C x C. We write its power series expansion in q x as 



r=0 

The H r {z) are holomorphic in z G C, the series converges normally on compacta of (x, G C x C, 
and Hq(z) = {q^ +z \ g 5 )^- Define new holomorphic functions by 

for r € Z> . Then fT («) = (<f +2z ; q 9 )^ and 

$r(x,z) = -5— ^ 5- — — > HJz)q rx 

with the series converging normally on compacta of (x, z) £ C x C. A direct computation now 
shows that 

= WJ{X ' Z ) y Hr (z)q rx . 

By (5.1) and (4.2) we furthermore have 

(<La;, " f ~ ?~* + f K + Q- 2K )$R(; z) = 0. 
Hence $r(x, z) satisfies the characterizing properties of 5>j(x, z). ■ 
By the proof of the theorem we in particular have 

/• f+2a; s\ 

(-? 5+ *;g 2 ) 00 

Substituting the explicit expressions of and ^fj in this formula gives the following quadratic 
transformation formula for very-well-poised 8^7 series. 
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Corollary 5.2. 

3 

(qfixz,-^-, , -g2axz;gj 



3 „ „2/W2 S W 7 \ -qzxz*;— , — , -az, qifiz, -q*x; q, -q^x 



[^x^-qixz^^-q)^ \ a P 

(- q ^f,-q^xz*,-<!f;qZ) 



\ a ' y /oo 

/ g/3xz 4 gV 2 2 2 2 gpz 2 gax\ 

x 8 Wt ; — 2~,-gz > _2: '9P z , ;<T, 5~ (5-2) 

\ a cr a p J 

if both \q*x\ < 1 and 1 < 1. 
Remark 5.3. 

(i) As observed by Mizan Rahman (private communication), the quadratic transformation 
formula (5.2) for $,Wj is equivalent to the quadratic transformation formula [6, (3.5.10)] 
by applying to the left hand side of (5.2) the transformation formula [6, (III. 23)] with 
parameters (a, 6, c, d, e, f) in [6, (III. 23)] specialized to 

( i o q*z i qz i \ 

I -q^xz , — —,-az,-q2x, —,Q 2 P Z I 

and to the right hand side of (5.2) the transformation formula [6, (III. 23)] with parameters 
(a, b, c, d, e, /, q) in [6, (III. 23)] specialized to 



qfixz^ _ 2 2 _ o2 __ 2 gV gpx 

v2 



-,-gz , -z ,g/3 z , 



a a* a 



(ii) The formal classical limit q ~\ 1 of (5.2) can be computed after replacing the parameters 
(x, z, a, f3) in (5.2) by (— x, q z ,q a , g' 3 ) and moving all the infinite g-shifted factorials in (5.2) 
to one side. The resulting classical limit turns out to be trivial, since both sides reduce to 

(l + 2z-2a } \ + 2z + 2f3 \x \ 
21 \ l + 4z ' {1 + xYJ ' 

where 2-F1 is Gauss' hypergeometric series (in contrast to the classical limits of the quadra- 
tic transformations of very- well-poised s^'s from [18, § 5], which reduce to nontriv- 
ial quadratic transformations for 2F1). On the polynomial level (5.2) reduces to the 
quadratic transformation formula for the continuous g-Jacobi polynomials, see Section 6, 
which is known to reduce to nontrivial quadratic transformations (see [1, (4.24)]) on the 
classical level when taking the limit q "f 1 after replacing the parameters (x, z, a, /3) by 
(-«*■ 5*, 5**, 0^)- 

There is a dual version of (5.2), which can be obtained by using the selfduality of <3?(x, z) to 
both sides of the quadratic transformation formula <3?j(x,z) = $r(x, z) before substituting the 
explicit expression as a series. It leads to the following quadratic transformation formula. 
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Corollary 5.4. 



(-q2axz,-q2p z , q -^ r ;q} ( q2 X 2 z qx 1 ax 1 q*z i 
a ~ — sW 7 77—; —7 ,-q'x,- — ,q'x, —;q,-q^z 

f_q*xht_. a \ V P ®P P P , 



(-qaPxz\^,-^f,^ q 2 )oQ 



(-q 3 x 2 z 2 ,-q 2 z 2 , q ^^;q 2 ) oc 

rxr ( 22 Q 2x QP X ax a IP 2 * 2 \ fK o\ 

x 8 W 7 -qx z ;—, , — —,qa/3x, ;q ,-qx . (5.3) 

\ ap a p a J 

Also (5.3) can be related to [6, (3.5.10)] by applying the transformation formula [6, (III. 23)] 
to both sides, cf. Remark 5.3(i). 



6 Polynomial reduction 

Both the asymptotically free eigenfunctions &(x, z) of the Askey- Wilson second order difference 
operator, and the Askey- Wilson function £(x,z), reduce to the Askey- Wilson polynomials 
when z is specialized appropriately (see, e.g., [10, 25]). Concretely, for 3>(x, z) we have for 

n G Z> , 

a~ 2n (a s b s , a s c s , a s d s : q s ) — l-;q s ) 
V ' (q(.»-Q>a a b a c a d a] q>) n St d (-K-v-ns) nK ' ' ' 

with P n (x) = P n (x; k, A, v, q, s) the normalized Askey-Wilson polynomial in q x + q~ x of deg- 
ree n, 

PM . , fq~ ns ,q {n - 1)s a s b s c s d s ,a s q x ,a s q- x s . 

\ &sUsi ^sC s , dgds 

This can be proved directly from the expression of <3?(x, — n — v — ns) as a very-well-poised 
sWj series (see Proposition 2.2) by first applying Watson's transformation [6, (III. 17)] with 
parameters (a,b,c,d,e,f,q) in [6, (III. 17)] specialized to 

/„2(1— n)s+x „s+x n (2—n)s „(1— n)s „(1— n)s „(l—n)s \ 

/ y y y y y y s \ 



a 



followed by Sear's transformation [6, (III. 16)] with parameters (a,b,c,d,e,f,q) in [6, (III. 16)] 
specialized to 

(q(l—n)s q(l—n)s q(l—n)s g2(l—n)s q(l—ri)s+x g{\—n)s—x \ 
a s d s ' b s d s ' c s d s ' a s b s c s d s d s d s ' J 

Formula (6.1) can also be proved by observing that P n {x) is an eigenfunction of T> x with eigen- 
value a(q n + 1) + a^ 1 (q^ n + 1) (see [1]) of the form 



{q- ns ,q {n - 1)s a s b s c s d s -,q s ) n (-a s rq^ n+ V 
p n{x) = — -p- 77; q nx + 2^ c m q mx , 



(a s b s ,a s c s ,a s d s ;q s ) n 
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hence up to a multiplicative constant it must be equal to <£(x, — k — v — ns) (this argument 
generalizes to the setting of Macdonald-Koornwinder polynomials, see [13, 25]). Also the Askey- 
Wilson function £(x, —k — v — ns) (n S Z>o) is a multiple of the Askey-Wilson polynomial, 

(a s b s , a s c s : q s ) 

E(x, -n-v-ns)= K 3 S J/^ y ; °° P w (x). (6.2) 

V a a d s ' * ) oo 

This can again be proved directly using transformation formulas and the expression of £ (x, z) 
as a series, see [10]. It can also be proved using the c-function expansion of £(x,z) (see 
Proposition 2A(iv)) and (6.1), since c(x; k + v + ns) = for n E Z>o- 

Specializing now z = —2k — ns (n G Z>o) in the identity < 5j(x, z) = z) and using (6.1) 

twice gives, after straightforward simplifications, the following result (we take s = 2 without 
loss of generality) , 

, fq- 2n ,a<f,aq-*,q 2n a 2 b 2 2 2 \ , fq- n ,aq*,aq-*,-q n ab 



ab,qa z ,qab ' ' J \ ab,q2a,-q2a 

for n € Z+. This quadratic transformation formula was first proved by Singh [22]. It was 
reobtained by Askey and Wilson [1, (4.22)] and interpreted as quadratic transformation formula 
[1, (4.22)] for the continuous q-Jacobi polynomial (to get it in the same form one has to apply 
Sear's transformation formula [6, (III. 15)]), see also [6, § 3.10] for further discussions. It was also 
obtained in [20, (3.20)] {—a 2 q in [20, (3.20)] should be —abq) by specialization of the quadratic 
transformation formula [20, (3.16)] for the R- function. A similar polynomial reduction can be 
done with the dual version (5.3) of the quadratic transformation formula, in which case it reduces 
to the quadratic transformation formula [1, (3.1)]. 

7 A theta function identity 

By the c-function expansion of the Askey- Wilson function (see Proposition 2A(iv)) and by 
Theorem 5.1 we have 

£j(x, z) = cj(x, z)<&j{x, z) + cj(x, -z)$j(x, -z) 
= cj(x, z)$ R (x, z) + cj(x, -z)$ R (x, -z) 

while on the other hand, 

£ R (x, z) = c R (x, z)$ R (x, z) + c R (x, -z)$ R (x, -z). 
The c-functions cj(x, z) and c R (x, z) are explicitly given by 

( „2k— z „— 2 _§+2A— z „f — re— A+z— 2. s\ 

yq > ~q i h z j ~q i > i ) 



cj(x,z) 



Wj(x, z)9(q- 2z , - q -2+*- x + x - q s) 



g(^-^)g(gI+A-f ; _ g I-*+*-|. g §) 

c R {x,z) = — — - — 3— j- -. (7.1) 

W R (x,z)0(q- z ,-qi +x ;q2) 

Note that Wj(x,z) = g (K+A+a>)(2K+*)/» = w R {x,z) and that 

cj{x,z) _ e(-q-2~ K ~ X+x - z ; q s )e{-gi +X --2,-ql +x - g§) 
cr(x, z) Q(-qi-z^ - q i ?a ^_ ? f-«+ x -§ . ? f ^ 

is not invariant under z — > —z for generic Hecke parameters, since it is zero at z = 2A + 
s + 2tt\/—1t but nonzero at z = —2X — § — 2tt^/—1t. Hence £j(x,z) cannot be of the form 



2 

a(x, z)£r(x, z) with a(x, z) meromorphic and s-invariant in the variable x. 
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Remark 7.1. Both £j(-,z) and £ R (-,z) are eigenfunctions of with eigenvalue g 2 + 

q~ z — q 2K — q~ 2K (for £ R (-,z) this follows from (4.1) and the fact that c R (-,z) is |-translation 

invariant). On the other hand, £ R (-,z) is an eigenfunction of with eigenvalue q^ + q~?, 

but this is not true for £j(-,z) since cj(-,z) is not ^-translation invariant. 

Remark 7.2. Ruijsenaars' i?-function lZ(-,z) [20] is a hyperbolic eigenfunction of the Askey- 
Wilson second order difference operator (implying in particular that it admits an analytic con- 
tinuation to the regime \q\ = 1) satisfying lZ(—x, z) = lZ(x, z). For the i?-function 1Z a quadratic 
transformation formula of the form TZj(x,z) = lZ R (x,z) holds true, see [20, (3.16)]. The dis- 
crepancy with the fact that £j(x,z) ^ £ R {x,z) is not unexpected since the i?-function 1Z in 
the trigonometric regime \q\ < 1 has a nontrivial factorization in Askey- Wilson functions: it 
expands as a sum of two terms, each term being essentially the product of two Askey-Wilson 
functions, see [3, Theorem 6.5]. 

The quadratic transformation formula &j(x, z) = & R (x,z) implies the following result for 
the connection coefficients in (2.8). 



Proposition 7.3. We have 

cj(x,-z) _ c R (x,-z) cj(x,z) - cj(z,x) _ c R (x,z) - 4(f,2x) 

4(z,-x) ~ 4(1,-2*)' c d (z,-x) 4(f, -2s) 

Proof. Using (2.8) and Theorem 5.1 we have 

' c,j{x,z) - c d j(z,x)\ cj(x,-z) 



(7.2) 



Cj(Z,—X) J Cj(Z,—X) 

c R (x,z) - c&(f,2x)\ c R (x,-z) 



V 4(1' ~ 2x > J c Uh~ 2x ) 

since both sides are equal to 3>j(— x, z) = & R (— x, z). By a straightforward asymptotic argument 
(compare with the proof of Corollary 2.8) it follows that the coefficients of $>j(x,z) (resp. 
of <3?j(x, —z)) on both sides of this equation should be the same. 
Alternatively, one verifies by a direct computation that 

cj{x-z)_ 2^ e(q 2x ,q^ +2X+z ,q 2K+z ;q s ) c R (x,-z) 



q 



X 



yielding the first equality of (7.2). Combined with (7.3) this implies the second equality 
of (7.2). ■ 

Since c(x, z) is s-translation invariant in both x and z, it follows from the right hand sides 
of (7.2) that the quotients 

cj(x,-z) cj(x,z) - c d j(z,x) 
c d (z, —x) ' c d {z, —x) 

are |-translation invariant in x (although cj(x,z) is not, cf. Remark 7.1). 
Using (7.1) and using the explicit expressions 

nj K+A— X K—X—X „f+K+A— x „% — k—\-\-z—x s\ 

c d l z % \ = ^ ' g 'g 2 » *? 2 >Q ) 

J ' W d (z,x)6(q- 2 *,-qi +z ;q s ) 



. , s 6*(fl K+A 2,-g K A 2 Q4 2,-04 K + 2 2-02 



x)6(q~ x , -g3+ A+2 ; gs ) 
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with Wj(z,x) = g («+A+a:)(2«+ z )/ s and W^(z,x) = g 2 ( K + 2 K K + A +! V s , we have explicitly, 

Cj(x,z)-cj(z,x) = 2x(2k+z)/s 0(q 2x ;q s ) 

0(g K + A+x _qK-X+x gf+K+A+x _^§-k-A+x+2. ^ 

J 0(g 2K " 2 , -g" 2 , g§+ 2A " 2 , -03-*-*+*-*, -gl + 2 ; g s ) 
X \ 0(g- 22 ,-gf +K - A+:c ;g s ) 

Q { „k+X— x _k— A— x _#+k+A— a: f — k— A— £+2. „s\ | 



-2x 



and 



•<;q s ) 



c R {x, z) - 4(|, 2x) _ q _ 2x(2K+z)/s 0(q 2x ; gi) 



4(|, -2x) 0( g «+A+s,-g«->+s, g f+ a: ,- g !-' t + a: +§ ;gf) 

'0(g 2K - 2 ;g s )0(gf +A -f,-gf- K+a: -f,-gi +A+ f;gf) 



X I 6»(g- 2 ,-g4 +:r ;gf) 

0(g K + A_a; , -g* - * -21 , gl _a; , -gf- K - a; +f ; gf ) 1 

^^f) J ' 

Hence (7.2) is equivalent to the following theta function identity. 
Corollary 7.4. 

2 t2 2 t 2 1 2 i2 j 2 2^ n fii 1^ i r 2 

a ,b c ,—b ,qb d , -, —go ; g — t) I o , -, aca, — —,qacd 



acd 1 ' J v ' ad' ' d ' 6 2 cd' 

0(-ga6 2 cd;g 2 )0(-o 2 ;g) / / 2 i g§o g§d \ 

' \ be, —be, a ,q?bd, , — ; g 

\ ac b 



0[qa 2 ; q 2 ^9[—q^abc; g) 

- 6 2 , -g5a, acd, -^,q*a, - V 1 ! (7-4) 



i 

ac i g2a 



d 6c 
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